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Abstract 

We have calculated the fermion contributions to the shift of the position of the 
poles of the massive gauge boson propagators at two-loop order in the Standard Model. 
Together with the bosonic contributions calculated previously the full two-loop correc- 
tions are available. This allows us to investigate the full correction in the relationship 
between MS and pole masses of the vector bosons Z and W. Two-loop renormaliza- 
tion and the corresponding renormalization group equations are discussed. Analytical 
results for the master-integrals appearing in the massless fermion contributions are 
given. A new approach of summing multiple binomial sums has been developed. 

1 Introduction 



In our previous paper referred to as I in the following, we presented the two-loop cal- 
culation of the bosonic contributions to the W- and Z-self-energies on the mass shell and 
discussed the main properties and possible applications. In the present paper II we extend 
this calculation to a full Standard Model (SM) calculation by including the missing fermion- 
loop and mixed contributions. Throughout the paper, we shall adopt the terminology calling 
"bosonic corrections" the one's represented by diagrams without any fermions and "fermionic 
corrections" the remaining one's given by diagrams exhibiting at least one fermion loop. 

In I, in particular, we have shown by an explicit calculation that up to two-loop order 
the purely bosonic corrections to the position of the pole of the gauge-boson propagators are 
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real for arbitrary values of the Higgs-boson mass. In contrast, since the W- and Z-bosons 
decay, at leading order, into light fermion pairs, the light (massless) fermion corrections 
give rise to a non-zero imaginary part at the one loop level already. The problem of gauge 
(in) dependence of the complex pole has been extensively discussed in literature 0. Only 
recently, two important results have been proven to all orders in perturbation theory: i) the 
position of the complex pole is a gauge independent quantity [Sj; ii) the branching ratios 
and partial widths associated with the pole residues are gauge independent jl]. Moreover, 
it has been shown that the pole mass of the W-boson is an infrared finite quantity with 
respect to massless photonic corrections. An alternative proof of the infrared finiteness of 
the two-loop bosonic contributions to the pole of the gauge bosons was presented in I by 
explicite calculation. It is based on the fact that within dimensional regularization [S], in 
d = 4 — 2e space-time dimensions, the singular 1/e terms, which regularize both ultraviolet 
(UV) and infrared (IR) singularities, are absent after UV renormalization of the position of 
the pole in the propagators. 

In the present paper, besides from completing our previous calculation by including the 
missing fermion contributions, we will discuss in some detail general features and technical 
problems which are specifically related to these contributions. 

The paper is organized as follows. In Section 2 we briefly reconsider the definition of 
the pole-mass of the massive gauge-bosons within the SM and remind the reader of some 
notation given in I. The required analytical results for the massless fermion two-loop master- 
integrals are presented in Section 3. In Section 4 we discuss the UV renormalization of the 
pole mass and the interrelation between our results and the one's familiar from the standard 
renormalization group approach. In particular, we performed several cross-checks of the 
singular l/e"^- and l/e-terms. General aspects as well as numerical results for the finite 
parts are discussed in Section 5. Some technical details and a number of our analytical 
results will be presented in Appendices. In Appendix A we present a set of non-standard 
binomial sums which are needed for the e-expansion of some of the hypergeometric functions 
entering the master-integrals. The one-loop fermion contribution to the pole masses of the 
gauge-bosons are reproduced in Appendix B. The reducible two-loop corrections obtained 
by mass renormalization of the one-loop fermion contributions may be found in Appendix C. 
In Appendix D the bare two-loop on-shell self-energy contribution of the massless fermions 
are given in exact analytical form. The corresponding contributions involving the top quark 
are presented in terms of the first few coefficients of the expansion in sin^ 9w and specific 
mass ratios in Appendix E. 

2 Pole mass: definition and calculation 

The position of the pole sp of the propagator of a massive gauge-boson in a quantum field 
theory is a solution for at which the inverse of the connected full propagator equals zero, 
i.e.. 
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where n(p^, ■ ■ ■) is the transversal part of the one-particle irreducible self-energy. The latter 
depends on all SM parameters but, in order to the keep notation simple, we have indicated 
explicitly only the dependence on the external momentum p and in some cases also m, where 
m is the mass of the particle under consideration. This can be either the bare mass or the 
renormalized mass defined in some particular renormalization scheme. 

Generally, the pole sp is located in the complex plane of and has a real and an 
imaginary part. By writing 

sp = M^ - iMT , (2.2) 

the real part defines M which we call the pole mass while the imaginary part is related to 
the width F of the particle. This is the natural generalization of the physical mass of a stable 
particle, which is defined by the mass of its asymptotic scattering state ^iTj. 

For the remainder of the paper we will adopt the following notation: capital M always 
denotes the pole mass; lower case m stands for the renormalized mass in the MS scheme, 
while mo denotes the bare mass. In addition we use e, g and gs to denote the U (l)em, SU{2)i^ 
and SU{3)c couplings of the SM in the MS scheme. 

In perturbation theory p.lj) is to be solved order by order. To two loops we have the 
solution^ 

sp = + ^W(m^ m^ ■ ■ ■) + m^ ■ ■ ■) + ^(^Hm^ m^, ■ ■ ■)^W'(m^ m^, ■ ■ ■), (2.3) 

which yields the pole mass and the width F at this order. II*-'^-' is the bare (m = mo) or 
MS -renormalized (m the MS -mass) L-loop contribution to 11, and the prime denotes the 
derivative with respect to p^. 

In this paper we show by explicit calculation at the two-loop level that the fermion 
contribution (including mixed terms) to the propagator pole sp of a gauge boson is a gauge 
invariant and infrared stable quantity. For more details concerning the tensor decomposition 
of propagators and abbreviations adopted we refer to Sec. 2 of I. 

In order to find the relationship between the poles of the gauge boson-propagators and 
the MS masses m|, m^ we have to compute the one- and two-loop self-energies for the Z- 
and VT-bosons at = m| and p"^ = m^, respectively^. 

For the calculation of the two-loop fermion contributions we again use the strategy de- 
scribed in our previous paper I. Let us present here some basic features. In order to be 
able to work with manifestly gauge independent parameter renormalization constants we 
have to include the Higgs tadpole diagrams. To keep control of gauge invariance we work 
in the gauge with three different gauge parameters ^w, iz and ^-y. However, in order 
to avoid additional mass parameters like y/^mw and y/^mz (ghost masses), we expand 

^Similarly, it is easy to find the on-shell wave- function renormalization constant Z2. Up to two loops it 
is given by the following equation [5]: 

= 1 - ^(l)'(m^ m2, . . .) - ^(2)'(m^ m2, • • •) - U^^\m\m^, ■ ■ " {m\ , • • •) , 

where is a bare or renormalized mass. For a recent discussion see also 

^There exist a number of programs for analytical and/or numerical calculations of two-loop self-energies. 
A selection may be found in .10) . 
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the original propagators in an asymptotic series at = 1. For the purpose of checking the 
gauge invariance of our results it turned out to be sufficient to keep the first two terms of 
the expansion. As an additional cross-check we actually kept one more term. 

For calculating diagrams with massless fermion loops we utilize Tarasov's recurrence 
relations [TT] (a detailed discussion we postpone to Sec. 3). For diagrams with top-quark 
and/or Higgs-boson propagators we apply asymptotic expansions with respect to the heavy 
masses. In these cases we firstly expand the propagator in the weak mixing parameter 
siia^ 9w = 1 — ^w/^'z ~ 0-25 and get rid in this way of mw (or rnz)- AH diagrams with 
top-quarks and/or Higgs-bosons are divided into several prototypes which are presented 
in Fig. ^ The large mass expansion has been performed with the help of the packages 
TLAMM ^21 or by the program described in ^31- In the heavy top limit, the leading 
and next-to-leading contributions to the gauge boson self-energies have been calculated 
in [T3l[T3] and [ini, respectively. 

As usual, in order to preserve gauge invariance of the regularized theory, we utilize 
dimensional regularization [S] for our calculation. In order to avoid spurious anomalies we 
adopt an anti-commuting 75 in dimensional regularization [SJ El and impose vector current 
conservation by hand in triangle fermion loops which contribute to the two-loop two-point 
functions ^H|- The axial vector current anomalies ^H] cancel for each fermion family by the 
lepton-quark duality |20|. In contrast to calculations which involve vertex- and/or box- type 
diagrams [211 > calculation of self-energies we do not encounter problems with violation 

of Ward identities (see also [221 ) • 

3 Master integrals for the massless fermion contribu- 
tions 

This part is devoted to the calculation of the two-loop massless fermion corrections to the 
pole masses of the gauge bosons. For this class of corrections it is possible to work out the 
exact analytical results without expansion^. 

In contrast to the previous calculations performed in [221, where results are expressed in 
terms of the non-minimal set of scalar two-loop integrals, we use here Tarasov's recurrence 
relations ^T] which allow us to reduce the number of integrals to a minimal set of master- 
integrals in our results. The analytical results for these diagrams are presented in [221 • Here 
we present an independent analytical calculation of the relevant master integrals shown in 
Fig. 121 Besides the known one's, we consider here two new master-integrals, shown in Fig. [3 
which contribute to the two-loop on-shell top-quark propagator in the limit of massless gauge- 
bosons. They may also be considered as the "naive" parts of the asymptotic expansion of the 
original diagrams with massive W- and Z-bosons in the limit, when the external momentum 
approaches the heavy mass-shell 

^In contrast to the bosonic corrections where aU diagrams can be expanded from the very beginning in the 
smah parameter sin^ 9w the individual diagrams with massless fermion-loops develop threshold singularities 
which behave like powers of Insin^ 9w- To control these terms we need the exact analytical result. 
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Figure 1: The prototype diagrams and their subgraphs contributing to the large mass 
expansion for two-loop diagrams with heavy propagators. Thick and thin lines correspond 
to heavy- and light-mass (massless) particle propagators, respectively. Dotted lines indicate 
the lines omitted in the subgraph. 



For the analytical calculation of these integrals up to the finite parts in e we use the 
method developed in [25] (see also (^Hl)- This approach is based on the possibility to retrieve 
the analytical results in terms of sums, as predicted by the differential equation method , 
from several of the first coefficients of the small momentum expansion |2HI- Whenever it 
is possible, we present the exact analytical results in terms of hypergeometric functions. 
Thereby, we apply Mellin-Barnes techniques developed in l3Uj. 

Let us introduce now the following notation for the finite sums which show up in the 
master integrals: 
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V2002 V1002 (ci,CF,,a,,|3) V2001 




FOOlOO VOOOl JOOO 

Figure 2: Diagrams corresponding to the master-integrals with massless fermion-loops con- 
tributing to the two-loop gauge-boson propagators. Bold, thin and dashed lines correspond 
to off-shell massive, on-shell massive and to massless propagators, respectively. 

Throughout the rest of this paper, the symbol Sa will always stand for Sa{n — 1) and Sa 
will stand for SaC^n — 1). The same notation applies for all the other type of sums, e.g., 
Va = Va{n — 1). We denote all master integrals as Tab-, where the first letter T = F,V, J 
indicates the topology in accordance with notation introduced in [TT]; indices A,B, - ■ ■ = 
0,1,2 characterize the relation of the corresponding internal mass to the external momentum: 
indicates a massless line, 1 corresponds to "internal mass equal to external momentum" and 
2 means that mass and momentum are different (see Fig. |21for details). In our normalization 
each loop is divided by (7r)^~'^r(l -|- e). 



6 



F20U0 



F00112 



Figure 3: Diagrams contributing to the two-loop on-shell top-quark propagator in the hmit of 
massless gauge-bosons. Bold, thin and dashed lines correspond to off-shell massive, on-shell 
massive and to massless propagators, respectively. 

• F00002 

For this integral we use the representation^ given in j2Sl 

■ i^oooo2 = 2C2 ln(l + z) + 2 \n{-z)U2i-z) + In'i-z) ln(l + z) + A ln(l + z)U2{z) 

-2Li3(-^) - 2Li3(^) + 2Si,2{z^) - 45i,2(^) - 45i,2(-^) , (3.4) 

where z = and Sn,p are Nielsen polylogarithms jHH IS2] • For the on-shell gauge-boson 
propagators two points in the variable z are interesting: (i) z = cos^ 9 w\p^ = m^,]; {u)z = 
1 / cos^ 6w [p'^ = rriz] ■ 

• F20200 

For this integral the following representation^^ is valid 

■ i^202oo = Lig (z) - 6C3 - C2 In y - ^ In^ y - 4 In ?/ Lia (y) 

o 

1^- ^ ,3 



+4Li3 (y) - 3Li3 i-y) + -Lh [-y') , (3.5) 



where we have introduced the new "conformal" variable 



Here we are interested in the case z = 1/ cos^ 6w only. 



^Another representation is given by Eq. (104) in j23j . 

^°The general case with two different internal masses is given by Eq. (106) in • 
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• ^2002 



The off-shell integral Vxooy with two different masses and unit powers of propagators was 
considered in [221 • The result, Eq. (46) of [221 '^^^ presented as a sum of F2 and F4 hyper- 
geometric functions of two variables. It turns out, however, that the extraction of the finite 
part from this result is difficult. The analytical result up to finite parts for this diagram (see 
Fig. [2} is given by^^ 



(M^)2^F2002 



1 1 1 

n 2 



1 



1 



2^2 2e 



-C2 + E^"- 



n=l 



11 . 1 

+ 4 

n 



2n 



1 1 



n 



2n 



-Wi + Si + 4 



1 1 

+ 



n + 1 n + 1 



V n+l J 



6Wi{n) + 4:Si{n) + 16 



(3.7) 



where z = p^/M'^ with M an internal mass. This result can be written also as 



p_ 



2\ 2e 



2002 



1 

'2? 



2 l-y 



19 

T 



rLi.f^) 



ln(l - z) 



1 + 1/ 



7 



Li2(-y) - -Li2( V) - 7 In' 2/ + 2 ln(l + y) In y - 4 In y + -(2 



(3. 



where y is defined in fj3.6|) . If we compare this with the cases considered previously in 
as a novelty, in Eq. ()3.7p we encounter series with shifted index of summation (n + 1 instead 
of n). The latter give rise to the appearance of terms of different weight in the finite part. 
They are connected with the presence of UV-divergent subgraphs. For this diagram only the 
point z = 1/ cos^ 6w [p"^ = is of interest here. 



• J 



002 



For the off-shell case, we have two master-integrals: one with all indices equal to one, and a 
second one with indices 1,1,2. Their e-expansion up to order 0(6:') are given by Eqs. (Al) 
and (A2) in Appendix A of |211. For our task only the following finite parts are needed: 



Joo2(l, 1, l,p',m') 
^002(1,1, 2,/, m') 



[m 



[m 



2\-2e 



1 z-Q ^ 13 1 

'2? + ^-^'-' + T^ + - 



1 1 



2z 



■In(l-z) -Li2 {z) 



2.= 2. ' 2 C,H-i^Mi-.)-Li,(.) 



For this diagram again two expansion points are important: z = cos' 6^ and z = 1/ cos' 6^. 



The general case with two different internal masses is given by Eq. (95) in |23|. 
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• Vioo2 

Let us present some useful details for this integral. Its Mellin-Barnes representation is 
^ioo2(a,(7i,(72,/3,p^,m^,M^^ ' - 



I 



T{l + e)T{l + e) 

r2(l + e)T{ai)T{a2)r{a)r{(3)T{d-ai-a2) 2m 
r(-s)r(/3 + s)r{a+l3+ai+a2+s-d)r{2d-2ai-2a2-2p-a-2s) 

r{f-ai-a2-a-p-s) 



(3.9) 



where P^'^\p, m) — 1 / {p^ — w? -\- \eY . Closing the integration contour to the left (s < 0) we 
find the following result for an arbitrary set of indices: 



\^-'^fV{a, + a2-mi-a,)Va-a2) 



X 



( _ ^ d-a-fS-ai -(T2 



m 

M2 



2 \ d—Oi—ai~a2 



T{P)T{a + ai + a2-d)T{2d-2ai-2a2-a) 



(i-(7i-cr2-f , d-ai-a2 + ^ 

— a — (Ti — (72, a — (Ti — (72 



Am 



M2 



r ( q; ) r ( q; + /9 + (7 1 + (72 — (i ) r ( (i — q; — (7i — (72 ) 

ni) 



r(f -q;-(Ti-(T2) 



1,^,a+f3+ai+a2-d 
l+a+ai+a2 — d, | 



X 



(3.10) 



We introduce the short notation 



1^1002(1, 1, 1, m^ M^) 1^,^^, = 



1002 



for the master integral a = P = ai = a2 = I- The result for this particular integral reads 



^1002 - 
2F1 



(M2) 



r(i-e)r(i + 2£) 



2(1 -£)(1 -2e)e2 r(l + e 



X 



2-e 



_ ^.-.2(l-.)(l-4e)r(l-4.)rtl-e) ^ 1. 1 - 2. 



(2 - 3£)(1 - Se) r(l - 3£)r(l - 2£) 



3-3e 



4xj 
(3.11) 



where x — Using the quadratic transformation for hypergeometric functions 



2-e 



4x = 



1 + VI - 4x 



1,£ 

2-£ 
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1 - y/l-4x^ 



(l + X) 2i^l 



1, e 
2-e 



X 



where x is defined in Eq. ()A.2|) . the first term can be reduced to a new 2-F1 function, whose 
all-order e-expansion is given by Eq. (2.14) of 



1, e 
2-e 



l-e 



2{1-2e)z 



1 + z- (1-z 



A~2e 



j=l k=l 



(3.12) 



Note, that up to order the expansion of the given hypergeometric function can be extracted 
from Eq. (A. 3) of ^j. The second term can be reduced to a 2-^1 function of the type 
considered in Appendix A (see Eq. ()A.22|) ). To this end the relation 



[1 + 2(a - b)e] z 
2(2 - he) 



3 
2 

Z-he 



1, 1 — ae 



;i - z) 



3 
2 

2 -he 



1, 1 — ae 



has to be applied. Combining the relations given above we may represent the integral as a 
series: 



(M2)2-ri002 



1 3 1°° f") 
2^2 2e e^^n+1 



+ E 



/2n\ 



-5Si{n) + ASi{2n) 



-\nx + 2Si{n) -2Si{2n) 



Inx + C2 



n + 1 



n + 1 



+3^2 (n) - 4^2 (2n) + 6Sl{n) - 10Si{n)Si{2n) + ASl{2n) 



^^S^jn) ^ S,i2n) ^ 4 



n + 1 ' n+1 ' {n + iy 
or in analytical form: 

.2 \ 2£ 1 



+ 4Si{n) - 4Si{2n) + 



n + 1 



(3.13) 



X 



1002 



2e^ 



+X In a; — In a; 4x — 1 + 



5 (1-X^) 
- + X In X ln(l + x) 

2 X 
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X 



[ln(l - x) + 21n(l + ;^)] - C2 - - + XC2 



X 



4 ln(l + x) - 2 ln(l + x) ln(l - x) - 3Li2 (-%) - Li2 (x) - 2 ln2(l + x) 



(3.14) 



The all order e-expansion for this diagram can be deduced for x > 1/4. Closing the contour 
of the integral ()3.9|) on the right (s > 0) we obtain an alternative representation. The latter 
also may be obtained from ()3.10|) by analytical continuation of the hypergeometric function 
to the inverse argument: 



-4e\l-2e){m') Vi 



2e 



'1002 



ly-^- r(l-e)r(l + 2e) / l,e 



Ax 
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2 r2(i -e)r(i + 2£)r(i -45) 



2^1 



;i - 3e) r(l + e)T{l - 2e)T{l - 3e) ^ ^ 
[1 - 2e) r3(l - 2s)r{l + e)r{l + 2e)ix [x 



1,-1 + 3e 
h + 2e 



Ax 



1 

1 

Ax 



(3.15) 



The e-expansion of the first hypergeometric function 2-^1 is given in |26j. It only contains 
log-sine integrals [32J. The second term at first may be transformed using the Kummer 
relation to a new function (see Eq. (4.13) in jSZj) whose e-expansion is given by Eq. (B.13) 
in |37^ and includes a Lsc function ()A.20|) . The analytical continuation of the e-expanded 
result to X > 1/4 can be constructed in the manner described in Sec. 2.2. of [H7] . 

The point of interest here is a; = cos^ 9w [p^ = ^w] ■ 
expansion parameter x will denote x = m'^/M'^. 



In the rest of this section the 



V 



2001 



For this integral we have the series 



1 3 1^ f") 



- In x + 25*1 (n) - 25*1 (2n) + 



1 



+ E 



2n 



x"< In X - 



3Si{n) - 2Si{2n) + 



n+1 



+ 2 



+2Sl{n) - 2Si{n)Si{2n) + 3 



Siin) ^Sii2n) 



n+1 n+1 (n + 1 



+ 



n+1 

Ins — 3V2{n) + 5*2 (n) 
- + ASi{n) - ASi{2n) + 



00 

n=l 



1 1 

+ 



1 



{n + ly n + 1 n 
and the analytical solution 

.2\ 2£ 



-C2- 



X 



Vn 



2001 



+X In X — liix 



1 1 
2^~s^ 



5 (1-X^) 

- + X In x ln(l + x) 

2 X 



4x + 



X 



ln(l-x' 



C2-y 



2X 

(1 - x') 

2x 



[41n(l-x) ln(l + x)-81n(l + x)]- 



'1-x) 



X 



ln(l — x) 



5Li2 (x) + 8Li2 i-x) - Li2 (x')J - — Li2 (x 



n + lj 
(3.16) 



(3.17) 



The physical point is x = cos^^vf- 



11 



• F20100 

The series and the analytical solution for this integral are given by 



M F20100 



and 



1 g f2n\ 



X 



+ 



n=l 



n n 



-- In^ X H — ln(— x) Inx + 
4 2 ^ ' 



2Si - 2^1 



n 



-Si + 5*1 

n 



In(-x) - C2 - - ^^2 + 3^2 



-352 + 65i5i-352-4^ + 4^--l-2- ^ ^ 



n 



n 2n^ 



'ft") 



(3.18) 



p'i^20ioo = 3 ln'(l + x) - 2C2 ln(l + x) - ln'(l + x) H'X) + In X H'X) ln(l + x) 



--ln2xln(l + x) + 21nxLi2 - 2Li3 (-x) - ^Lis (x'j + -Lis (x) - ttLis (x 



6 



(3.19) 



respectively. The physical point is again x = cos^Ow- 
• ON-SHELL2 

We are interested also in the integrals -foo2oo(l, 1, 1, 1, 1,^^,"^^), Vooo2(l, 1, 1, l,p'^,m'^) and 
Jooo(l; 1; 1;P^) when the external momentum is on-shell. Up to finite parts these diagrams 
have been calculated in for arbitrary values of the external momentum. For our purpose 
it is important, that the limit p"^ — exists and is smooth, such that we can substitute 
the on-shell values of these diagrams^^ without problems: 



Fooioo(l, 1, 1, 1, 1, m , m ) = — 



Voooi(l, 1, 1, l,m^m2) = {m^y 
Jooi(l, 1,1,^2) = (m2)i-2- 



-3C3 + 2mC2 



-2e 


1 


5 




~2^ 


~ 2e 


1 


5 


11 


2? 


~ 4e 





19 

Y 



2C2 



Jnnn(l, 1, l,m 



2\l-2e 



1 13 1 

1 iir- 

Ae 18 2 



^^We note that the sign of the imaginary parts of the foUowing diagrams: 
(Fooioo, Voiio, Vooio, Voooo, Jooo) collected in Appendix C of must be changed in order to get 
the correct answer. 
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• F20110 



This master-integral shows up as a Higgs contribution to the pole of the top-quark propagator 
in the approximation of massless W- and Z-bosons and a massive Higgs-boson {m'jj = 
M^, = ^ M"^, ^D- It may be obtained also as the leading term of a Taylor 
expansion of the original diagram with massive W- and Z-bosons in the limit, when the 
external momentum is on the heavy mass-shell For this integral we find the following 
series representation 



M F20110 



1 g f2n\ 



X 



n=l 



n n 



- In^ X Inx Co 

4 n 2^ 



I _ ^ - -2 3 1 

'2^'^ ~^ S2 — Si + 25'iS'i — 5;^ + - — 



(3.20) 



and the analytical result 



m'F2oiio = -C2 ln(l + x) + 2 X ln(l + x) + 2 In xU^ (-x) - SLig (-x) • (3.21) 



The on-shell result coincides with 



00112 



This master-integral again is important for the calculation of the Higgs correction to the 
pole-mass of the top-quark. Its representation is 



M^Fooiia = -El 1— <!-ln'x + 



-2Si + 2^1-3 



1 



n 



Si „ 5*1 



+^2 - 2^2 + 2St - ^SiSi + 2St + 6— - 6— + 6— 



n 



n 



(3.22) 



or, in analytical form. 



m'Fooii2 = 2C2 ln(l + x) + In' X ln(l + x) + 2 In xLi2 (-x) 



(3.23) 



The on-shell result again coincides with the one given in [HH] . 



4 Renormalization 



The pole mass is a gauge independent and infrared stable quantity. For the SM this has 
been shown to all orders of perturbation theory jS]. In order to calculate the pole mass 
in terms of renormalized quantities at the two-loop level we need to calculate the one-loop 
renormalization constants for all physical parameters (charge and masses) as well as the 
two- loop mass renormalization constant. Not needed are the wave-function renormalizations 
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and the unphysical sector renormalizations. In order to obtain a gauge invariant result in 
the SM, however, we have to add in a proper way the tadpole contributions 1^0] . 

We first perform the UV-renormalization within the MS scheme in order to obtain finite 
results. In a next step we work out the relation between the on-shell and MS parameters. We 
adopt the convention that the MS parameters are defined by multiplying each L-loop integral 
by the factor (exp(7)/47r)^^. As usual we denote by fi the MS renormalization scale. For 
the RG functions we use the following definitions : for all dimensionless coupling constants, 
like g, g', gs, e. A, the /3-function is given by p-'^^g = /3g and for all mass parameters (a mass 
or the Higgs v.e.v. v) the anomalous dimension 7^2 is given by /i^^lnm^ = 7^2. 



4.1 One- loop charge renormalization 

The relationship between the bare charge Cq and the MS charge^^ e reads 



with 



Co 



7(1) 



167r2 



7 2 

1 — rip 

2 3 



9 ' 



(4.24) 



(4.25) 



where np is the number of fermion families and Nc the number of colors. In the SM we have 
Up = 'i and Nc = 3. The corresponding /^-function reads 



_d_ 



.3 , 



167r2 



7 2 

1 Up 

2 3 



9 



+ 0(e^ 



(4.26) 



where g and g' denote the SU{2)l and f/(l)y gauge couplings, respectively. 



4.2 Mass renormalization 

The mass renormalization constants {Zm)v at two loops may be written in the form 

my^, - myif,) {Z^)y - myif^) \1 + ^^Zy +^^^^Zy +j^^^,Zy j, 

(4.27) 

where V stands for any of the bosons Z, W or H. We shall use the same notation for 
fermions V ^ f, where we are interested in particular in the top quark f = t. In addition 
to the masses, we have one coupling constant as a free parameter of the SM which we have 
chosen above to be the electric charge strength e = gsinOw The one-loop mass counter- 
terms are well known [10] • We divide all corrections into bosonic and fermionic parts, 
Zy = ^^/boson -^yfermion- '^^^ purely bosonic contributions have been given in I. The 
fermion-loop corrections yield the following MS renormalization constants 

^•^AU MS-parameters, like e, g' are ^-dependent quantities. 
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(4.30) 
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(4.31) 



where sums run over all leptons (e,yU,r) (leptons) or quarks {u,d, s,c,b,t) (quarks). All 



masses here are MS masses and depend on the renormalization scale fi: 



m 



V 



m^(/x). 



In general, we will divide the two-loop fermion contribution to the mass-renormalization 
constants into two parts: the first one includes the contribution from A^^ massless fermion- 
families and the second one the contribution of Nh lepton-quark families exhibiting one 
massive u-quark of mass rrit (massive top-quark). The number of fermion families is n^? = 
+ Nfi, where in the SM we have = 2 and = 1. The results are given by the 
following exact expressions 
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^(2,1) 
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The renormalization group equation may be utilized to verify the higher order pole terms 
l/e^. It is given by Eq. (4.16) of I 



The value of [3^p = ("ft + 1)^f) may be calculated from the relation 



/3™ =in9„, = ip^^ f (ZS'" - 4"') + , 
" 2 smfc'iy yloTT^y ^ ^ 

where Z?^^-' is given in ()4.26p . and cos^^vf = ""^vv/^l- additional relation which holds for 
the 1/e-terms is the following (see Eq. (4.19) in I): 



sin* ^w- ^ sin^ ^vy cos^ ^vk = - ( tIz^ ) ( ^i^'^ - ^J'^ ) cos^ dw • 



The two-loop /3-functions for g and are given in and read 

1 10 \ g'^ (\ 95 \ g'^ /3 1 \ c/'V 



1^-3 = I + Y-i^j + U + 54"^; (16^ + U + 2"^; (16;^ 

17 ^('^^(^ 
'24m2^(167r2)2' 

. I , 43 2 \ / 259 49 ^ , ^1 , 



12 3 ) 167r2 V 12 6 7 (167r2)2 V4 6 / (167r2)2 

2 

t 



3 m2 



2^2 



8 (167r 

In all cases we could verify our results to satisfy the above RG equations. We also have 
established that our two-loop RG equations ()4. 32^ - 1)4. 39p . calculated in the broken phase, 
are related to the ones found in the unbroken theory a long time ago On the one 

hand the relationship between the UV counterterms and the RG equations in the broken 

^^Only recently, these results have been confirmed in 143| by an independent calculation. 
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and the unbroken phases for massless couphng constants is a trivial consequence of the Higgs 
mechanism (spontaneous symmetry breaking), which by definition must preserve the Ward- 
Takhhashi- and the Slavonv-Taylor-identities. On the other hand the Higgs field vacuum 
expectation value v appears only in the broken phase and its renormalization is looking very 
different from the ones one has to perform in the symmetric phase. In the broken phase 
the MS parameter may be defined by the ratio f^(yU^) = 4m^(/i^) sin^ 6'vi/(yU^)/e^(/i^) 

and the MS Fermi constant^^ by \/2Gp{ii^) = The relationship with the unbroken 

phase follows if we define, alternatively, v'^{iJ?) = rn? {iJ?) / , where rn? and A are the MS 
parameters of the scalar potential in the symmetric phase. As a consequence the following 
relations hold between the RG functions of the broken and the unbroken phase: 

Iv^ =lw - {iw - Iz) - 2— = 7^2 - ^ . (4.40) 

sm d e A 

Our results allow us to confirm the /3-function for the combination 7^2 — (3\/\, not however, 
the two-loop RG equations for A and im? independently. Given the RG equation for f^, we 
are able to write the two-loop RG equation for the top-quark and Higgs masses using the 
two-loop result |121 for the corresponding coupling, i.e. 7^2 = 7^2 (for details we refer 
to gg). 

These considerations also shed light on the role of the tadpole contributions in param- 
eter renormalizations. The zero momentum transfer Higgs propagator which multiplies the 
tadpole loops implies (in the perturbation expansion) non-analytic terms proportional to 
l/rrijf ~ 1/A. In view of (j4.4Uj) and the mass coupling relations, like mw = gv/2 etc., it is 
obvious that such terms must appear in mass counter-terms (see above and Sec. 4.2 of I) as 
well as in the coefficients of the corresponding /5-functions. 



4.3 QCD corrections 

At the two-loop level there are also diagrams exhibiting gluon exchange which are contribut- 
ing to the gauge-boson propagators. The exact analytical result for the 0{aas) contribution 
of quark-pairs to the vacuum polarization function of the gauge-bosons have been calculated 
in |46j (see also [HI). The corresponding set of IPI diagrams yields a gauge independent 
and infrared stable (with respect to taking the limit of vanishing fermion's masses) contri- 
bution^^. For the sake of RG invariance, in contrast to the usual conventions, we include 
tadpole contributions which cancel in physical observables like Ar, Ap etc. jlH] (see also [T3]). 
For each quark species, there is one 0{aas) tadpole diagram which is gauge invariant. For 
a quark with mass rriq its bare contribution to the location of the pole of the gauge boson 

^^Note that our definition of the MS Fermi constant is different from the one used in |44j . 

^^It is interesting to note, that in contrast to the two-loop electroweak massless fermion contributions, 
where it suffices to expand all master- integrals up to the finite parts, the master-integrals J002 and J022, 
showing up as additional basic integrals when considering the QCD corrections to the W- and Z-propagators, 
in this case must be expanded up to terms proportional to and e, respectively. The situation here is similar 
the one discussed in |34| . 
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propagators reads 
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with Cj = 4/3 in the SM. Taking into account the 0{as) term of the MS one-loop top-quark 
mass renormahzation constant 
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(4.41) 



where = 9^/4:71 and the first five quarks are treated as massless. The terms proportional 
to mf come from the tadpole contribution and will cancel in observable quantities. Again, 
the renormahzation group equations can be used for cross-checking the l/e^- and 1/e-terms. 
The l/e-part satisfies the following relation (see Eq. (4.19) of I): 
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and the -Z'^^j denote the l/e-' -parts of the Zy"" renormahzation constants ()4.41|) . The terms 



proportional to 1/e may be calculated from the relation 



— 1 1 m: 
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(4.42) 



where the definition of Zy'^'^ has been given in ()4.27|1 . Here, we have to take into account 
that, in the presence of several dimensionless coupling constants, Eqs. (4.14) and (4.15) in I 
have to be modified accordingly: 
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Note that our Eq. (j4.42|) explicitly reveals, that the systematic inclusion of the tadpole 
contributions is important for the self-consistency of the renormalization group equations. 



4.4 MS renormalization of the propagator pole 

The calculation of the one-loop MS renormalized on-shell amplitude fly^ is well known (see 
e.g. |in|)- We get it by rewriting the bare expression in terms of MS parameters 
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(4.45) 



As explained in Sec. 5 of I, we may avoid the consideration of wave-function renormalization 
as well as the renormalization of the ghost sector and of the gauge parameters if we look 
directly at the full two-loop MS renormalized on-shell amplitude^''. The latter can be written 
in the form 
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where the sum runs over all species of particles j = Z, W, H, t and 
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The functions Zy'^\Xy^ and X^^y are defined in 

Appendix C of I, respectively. The derivatives of x'^y with respect to the masses may be 
found in Appendix O 

It is interesting to observe that (j4.45p and (j4.46j) account for the renormalization of the 
real and the imaginary parts, i.e., for the mass and the width, simultaneously in a unified 
form. 
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5 Results and discussions 

We have calculated the location of the pole in terms of the MS -mass for the massive gauge- 
bosons Z and W . As in I we will write the result in the form 



where both e and sin6'vy are to be taken in the MS scheme. The one-loop coefficients Xy 
for Z, W and H are known of course as exact results. We wrote them down for completeness 
in Appendix B. The light (massless) fermion contributions to the coefficients Xy are given 
in Appendix O in exact analytical from in terms of the master-integrals discussed in Sec. El 
Alternatively, these contributions may be represented in a sin^ ^v^-expanded form^^ in the 
same way as the other contributions. In the latter case, after summing all diagrams with 
massless fermion loops, all singularities of type In-' sin^ dy/ cancel which infers the infrared 
finiteness of the massless fermion contribution to the pole-mass. When the top quark is 
involved we perform an expansion and present the coefficients as series with respect to three 
small parameters. One "small" parameter is the weak mixing parameter sin^ Qy/-, the others 
are the mass ratios niy/m'jj and my/rn^ where V = W ot Z. The analytical values of the 
first three coefficients are presented in Appendix E. 

Again we have checked, that the relations between the pole- and the MS -masses, when 
calculated in the so-called modified MS scheme (MMS), coincide with the ones obtained in 
the ordinary MS scheme. 

We have developed and implemented into a computer program an efficient algorithm for 
the calculation of the massive Feynman diagrams with external momentum on the mass 
shell of one of the internal particles. The necessary analytical continuations which allow 
us to cover the hole region of parameter space has also been constructed and implemented 
numerically. 

After UV renormalization the pole-mass is a finite expression in the limit e —>■ 0. Since 
also the IR singularities have been regularized by dimensional regularization, this result 
implies the infrared finiteness of the SM contributions to the pole-mass. As a consequence, 
the pole-masses of the gauge-bosons are infrared finite quantities. Note that, in order to 
establish the gauge invariance of the location of the pole Sp, the tadpole contributions have 
to be taken into account. Also RG- invariance requires the inclusion of the tadpoles. 

Our calculation proves that the MS renormalization scheme, comprised in ()4.46|) . is self 
consistent and works properly in case of unstable particles. 

In contrast to the mass of a stable particle, the definition of a mass of an unstable particle 
is not unambiguous. The natural definition by y/Resp (see 12.11^?^ . which we adopted in 
this paper, is not always used. In particular LEP/SLC experiments have adopted another 
definition, which has been used by most of the theoretical papers on the Z line-shape before 
LEP/SLC started to take data (see for a pre-LEP status report). As a consequence 
the Z and TV-masses determined by LEP/SLC experiments and listed in the particle data 

^^Details concerning this expansion can be found in j45|. 




(5.1) 
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tables correspond to My = My + Ty rather than to the pole-mass My In this case the 
width r'^ entering the Z line-shape must be defined by linllz{s) ^ sTz/Mz and hence is a 
function of the cm. energy square s. While at one-loop order the width is determined solely 
by the gauge invariant fermion contributions at higher orders a s-dependent widths, which 
is not an on-shell quantity, will cause troubles with gauge invariance. However, one still may 
understand the LEP measurements to be defined from the quantities defined in terms of the 
gauge invariant position of the pole sp by the gauge invariant relationship My = My + Ty. 
In the latter case the mass may be considered to be defined via My /My = \sp \ /Res p. 

Our results for the two-loop mass renormalization constants in the on-shell and the MS 
scheme can be applied for the calculation of physical quantities in both of these schemes at 
the two-loop level. Examples of such calculations, where results of our paper I have been 
used, are the computations of the bosonic two-loop contributions to the muon life-time 
(usually encoded in the correction Ar [SI]) presented in [S^ . 

The numerical evaluation of our results is of interest in several respects. First of all, even 
though the calculated quantities are not observables, but only a subset of contributions to 
those, their numerical magnitude is interesting. It certainly sheds light on the size of possible 
corrections as well as on the convergence of the perturbation expansion. Last but not least, 
the numerical evaluation by itself is a challenge: the amplitudes are of remarkable complexity 
and the numerical stability and the efficiency of the evaluations are far from trivial. Although 
we are largely working with series expansions already, and not with numerical integration 
of the master integrals, we have to apply numerical multiple precision techniques to be sure 
that we get the correct answers. To this end numerical evaluations at the present stage have 
been performed in MAPLE were we can work at the desired precision only on the expense 
of CPU time. 

Since we are calculating the on-shell counter terms for the gauge bosons, it is natural 
to consider the on-shell scheme with the fine structure constant a and the masses as input 
parameters. For our numerical calculation we have chosen the following input parameters: 
a = 1/137.036, Mw = 80.419 GeV, Mz = 91.188 GeV. Neutrinos, leptons and quarks are 
taken in the massless approximation and the top-quark mass is taken to be = 174.3GeV. 
For the strong coupling constant we are using a fixed value = as{Mz) = 0.1185. Often 
the two-loop corrections are of comparable size to the mass effects of the "light" fermions 
at one-loop. We then take the following fermion masses for the numerical evaluation: 



Fermion masses in GeV 



me 


= 0.00051 


rriu 


= 0.003 


rrid 


= 0.006 




= 0.1056583 


rric 


= 1.2 


rris 


= 0.120 


rrir 


= 1.77703 




= 174.3 


rrib 


= 4.4 



The general features of our expansion have been discussed in I. While the series expansion 
in sin^ 9w for the known fixed value converges very well, the expansions in My/m^j for large 
iriH leads to a strong coupling problem and perturbation theory becomes useless above 

^^Note that M'^ — Mz — 35MeV and thus in order to obtain the pole-mass from the "experimental" 
Z-mass one has to subtract about 35 MeV. 
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about 800 GeV. In general, as is illustrated by Figs. El and El the first three coefficients 
yield a good approximation as long as the asymptotic series behaves well. Of course, the 
expansion in My/m'jj for smaller values of rriH starts to diverge below about 120 GeV. We 
have calculated six terms in each of the different expansion parameters, but only write down 
three of them in Appendix E Figs. EHHl show the corrections Aw = M^/m^(Mvi/) — 1 and 
Az = M'^/m\{Mz) — 1, respectively, as a function of the Higgs mass for intermediate 
and for heavy Higgs masses. The light fermion contributions are small relative to the other 
corrections. The same is true for the QCD corrections for Higgs masses above about 300 
GeV. In Fig. [TUlthe imaginary part and the absolute value of the pole position sp are shown. 

Very often the inverse of ()5.1|) is required. To that end we have to solve the real part of 
(12. 3|) iteratively for my and to express all MS parameters in terms of on-shell ones. 

The solution to two loops reads 
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where the sum runs over all species of particles j = Z, W, H, t and 
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I in the MS scheme and parameters 



stands for the self-energy of the jth particle at = m 

replaced by the on-shell ones. The analytical results for the derivative term J2 
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can 



be extracted from the results of Appendix C. In ()5.2|) my — My 



(5My)gj^j^g represents the 



finite part of the on-shell mass counter-term. The corresponding bare one is obtained by 



5my 



m^y - m 



adding the MS mass counter-term (see section 4.2' 

The relation ()5.2p involves a change from the MS to the on-shell (OS) scheme also for 
the electric charge. Let us consider therefore, in details, the relationship between the fine 
structure constants and aos = 

For the electroweak couplings we have to calculate the MS versions from their com- 
monly used on-shell values^^. The MS version of the fine structure constant is defined as 
a solution of the renormalization group equation (see I4.26|) and can be calculated from the 
UV-counterterms of the electrical charge e (see (|4.25p ). In perturbation theory the naive re- 
lation between MS and on-shell values of a is determined by the Thomson limit of Compton 
scattering and can be written as 
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(5.3) 



^'^The full set of coefficients may be found at jhttp : //www- zeuthen. desy . de/| ''kalmykov/pole/pole .html 

^^The QCD coupling is parameterized almost always in the MS scheme. Therefore ^^(/x^) usually is 
directly determined experimentally by fitting data to a suitable perturbative QCD prediction in terms of 
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where the sum goes over all fermions / with N^j = 1 for leptons and N^f = 3 for quarks 
and Qf is the fermion charge equal to 1 for leptons, -1/3 and 2/3 for d- and u-quarks, 
corresp ondingly. 

The terms proportional to logarithms of the fermion masses originate from evaluat- 
ing the finite part of the derivative of the photon propagator at zero momentum transfer, 
which we denote as 11:^^(0)^^. However, the low energy contribution of the five light quarks 
{u, d, s, c, b) cannot be calculated in perturbation theory. The free quark loops are strongly 
(non-perturbatively) modified by the strong interaction at low energies. In order to evaluate 
11;^^ (0) we may write it as 



n/^(0) = Ren^,(g2) - 
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Ren^^(g2)_n/,(0) 



77^ 



(5.4) 



where is chosen sufficiently large (typically M|) such that the first term on the r.h.s. can 



be calculated in perturbative QCD. In the limit of large momentum transfer s ^ m'l the 



result is 



Ren;J^(s) 



(5.5) 



M|) 



which we use in the relation ()5.3|) for the perturbative light quarks contributions (s 

On the grounds of analyticity and unitarity, the second non-perturbative term in ()5.4|) 
can be determined by evaluating a dispersion integral over the known experimental e^e~ — >■ 
hadrons cross-section (for details see 1^2] )• Usually, the cross-sections are represented in 
terms of the cross-section ratio 



R{s) 



r 



hadrons) 



where crfe^e 



r 



at tree level. In terms of R{s) we obtain 



Ren^^(g2) - n^T,(o) 



4g2 
~3 



/ ds 



R{s) 



s [s — q 
Ana X 



The shift in the fine structure constant in then Aa{s) = Ana x j^AIl^^{s). Using the 
experimental data for R{s) up to ^/s = Ecut = 5 GeV and for the T resonances region 
between 9.6 and 13 GeV and perturbative QCD from 5.0 to 9.6 GeV and for the high energy 
tail above 13 GeV one gets 



^"hadrons 



0.027572 ± 0.000359 



a 



-\M: 



128.952 ± 0.049 



(5.6) 



at Mz = 91.19 GeV. For numerical estimations we use the results of [SI]. 

The shift in a calculated by the dispersion relation corresponds to the on-shell scheme. 
Accordingly, since 11*^^'' depends on e by an overall factor only, we have 
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At zero momentum fermions contribute to the charge only via n^-y(O) 
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9e 

where 



and 



5AaMscb(M|) = -— In. 



977 V""" /i2 3 J 

the perturbative subtraction-term (j5.5p form the 5 hght quarks. The quark loops contributing 
to the on-shell gauge boson self-energies are evaluated here at a high energy scale and hence 
it should be save to calculate them in perturbation theory^^. 

Finally, we analyze the Higgs mass dependence of the MS parameter sin^ 9 defined as 

9 + 9^ 

The Fig. ITT] shows the correction to Sg^j^2 q = sin^ 9^p/ sin^ 9^ — 1 as a function of the Higgs 
mass. Since the unphysical mjj terms drop out by virtue of Veltman's screening theorem, 
the corrections do not blow up so dramatically in the strong coupling regime of large Higgs 
masses. In the region displayed between 150 and 800 GeV the two-loop correction stays 
below about 0.6% and has a sign opposite to the one-loop correction. The latter is one 
order of magnitude larger. 



6 Conclusion 

We have calculated the full SM two-loop radiative corrections to the pole masses of the 
gauge bosons W and Z. A number of conceptual problems dealing with the renormalization 
of unstable particles could be studied by explicite calculations: i) The position of the complex 
pole sp of the gauge boson propagators is manifestly gauge invariant after taking into account 
the Higgs tadpole contributions, ii) The renormalized on-shell self-energies are infrared 
finite, ill) Our calculation proves that the MS renormalization scheme, comprised in ()4.4(ij] . 
is self consistent and works properly in case of unstable particles, iv) Up to two-loops 
we explicitely confirm that the UV singularities and the related RG equations of the broken 
phase are completely determined by the unbroken phase i45 . v) The inclusion of the tadpoles 
is also required from the point of view of the renormalization group invariance. vi) Our 
results for the 2-loop mass renormalization constants in the on-shell and the MS scheme can 
be applied in calculations of physical quantities in both of these schemes at the 2-loop level. 

^■^Possibilities to evaluate them by non-perturbative methods via dispersion relations have been discussed 
in [HS]. 
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vi) A compact expression for the massless fermion contributions in terms of an irreducible 
set of master-integrals in given, vii) A new technique for calculating the e-expansion of 
some types of hypergeometric functions has been developed. 

Note that all results concerning general properties (like gauge invariance etc.) have been 
checked analytically. Since exact analytical results for many of the multi-scale master- 
integrals are not yet available, in this work we had to resort to asymptotic expansion tech- 
niques. The exact results in terms of a basis of master-integrals will by given elsewhere. For 
the latter one dimensional integral representations are available which are suitable for the 
numerical evaluations ^U]- However, for numerical evaluations in the heavy Higgs region 
above about 200 GeV, the computation in terms of our series expansions is numerically much 
more stable and much more efficient. 
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Appendix 

A Sums with binomial coefficients in the numerator 

Let us present in the following some useful formulae for the binomial sums which occur in the 
master- integrals ()3.18|) and ()3.20p . The sums which involve only harmonic coefficients, were 
investigated in [^2], while the sums with binomial coefficients in the denominator^^ were 
considered in fl^ IHUj . In this Appendix we present some results for sums with binomial 
coefficients in the numerator. A sum Fa of weight a has form 



fc=i 



where f{k) is a product of some finite sums (e.g. Sa, Sb, etc.). It is convenient to introduce 
a new variable x related to x via 



Y = , x=- —, V 1 — 4x = . A. 2) 

A general relation for sums of type ()A.1|) is then given by 

Fa-i{x) = x^Fa{x) = \±^x^Fa{x)- (A.3) 
dx l-x dx 



^"^In sums of this type were called multiple binomial sums. Here we will use the more adequate name 

inverse multiple binomial sums. 
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Relevant analytical results for sums of the type considered are: 



E 

k=l 



E 

k=l 



E 

k=l 
oo 

E 

k=l 
oo 

E 

fe=i 

oo 

E 

fc=l 

oo 

E 

k=l 

oo 

E 



k=l 



k / k^ 



'2k\ 
k I P 



-8Si,2(-x) ln(l + x) + 4Li3 (-%) ln(l + x) - 4Li2 (-x) ln'(l 

ln^(l + x) - Li4 (-x) + 4S2,2(-X) - 8Si,3(-x) , 
Si = ~ ln=^(l + x) + ^Li3 (x') - 2 ln(l + x)Li2 (x' 
+4Si,2(x)-2Si,2(x' 



2k\ - 2 

J ^^1 = -3 + X) + 2Li3 (x) - 4 ln(l + x)Li2 (x) 



'2k\ x'' 
k T^' 



+4Si,2(x) + 6Si,2(-x)-2Si,2(x'), 
= ln3(l + x) - 41n(l + x)Li2 (-x) " 4Si,2(-x) 



2k\ X 



—St = - ln^(l + X) + 8 ln(l + x)Li2 (x) + 4 ln(l + x)Li2 (-x) 



k 



-4Si,2(-x)+4Si,2(x') 
2 



'2k\ _ 

^ j -S,S, = - ln\l + x) + 6 ln(l + x)Li2 (x) - 6Si,2(-x) 



'2k\ x^ 
k ] T 



+4Si,2(x) + 2Si,2(x' 
2 



Z ln3(i + ^) _ 41n(l + x)Li2 (x) - 881,2 (x) 



2k\ X 
k j ~k 



S2 — S^ 

= -Li3 (x) - 2Li3 (-x) + ^Li3 (x') - Li3 (x) . 



All these low- weights sums may be obtained by applying ()A.3|) . 

The following representation is often useful for sums of the type just considered: 



X) 

(A.4) 

(A.5) 

(A.6) 
(A.7) 



(A.8) 

(A.9) 
(A.IO) 

(A.ll) 



E O '^m 



n=l 



1 



n / n 



(a -2)! Jo 



n=l 



n / n 



In X — In z 



a-2 



dz 

z 



(A. 12) 



where f{n) is an arbitrary combination of finite sums. For particular cases, f{n) = 1, 5*1, 6*1, 
S2, SiSi, S2 — Sf or V2, the corresponding expressions extracted from ()A.4|) - ()A.11|) allow 
us to perform the sum under the integral in ()A.12|) exactly and to obtain one-fold integral 
representations for a large class of binomial sums. For < x < 1/4 , a further simplification 
is possible by substituting, z = ^ sin^ 0. It leads to the following representation: 



^f2n\x- 



n=l 



(a-2)! Jo 



arcsin v 4a; 



COS< 
) 

sin ( 



ln(4x) — ln(sin^ 



a-2 



X 
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E 



C) 


fin) 


1 sin^ (j) 






n 


V 4 





= arcsin -\/4x) , (A. 13) 

where T,{f,6) is the short notation for this sum. Unfortunately, we could not find corre- 
sponding analytical results for these sums for arbitrary x {x < 1/4) in terms of known 
functions, like generalized Nielsen polylogarithms. Not even for the simplest case f{n) = 1, 
which yields 



n=l 



la+l 



arcsin V4x I 

In cos - 
\ 2. 



ln(2v^) - In sm I 



a- 2 



d In sin 6 



a solution could be found. The problem here may be compared with the one encountered in 
the context of the inverse binomial sums, considered in |56[ 15 7j. 
For X > 1/4 the representation 



E 

n=l 



n J rr 



TT 



+ 



'a -2)! 



TT 



4 n=l 



2 arccos 



2n\ 



n / n 



fin) 



In X — In 2; 



(a -2)! Jo 

00 

E 



cos 



2 L 



a-2 , 

z 

In + In ( 2 cos - 



n a-2 



X 



n=l 



2n^ /(n) 
n 



n 



4 cos^ 



(A.14) 



is available, where we have introduced a new variable z = l/(4cos^|) and S (/if) is a 
constant, related to some (combination of) hypergeometric functions pFq of argument unity. 
For each given choice of f{n) these constants are different, but in all cases they are expressible 
in terms of "even basis" elements (see jSH] and Appendix B.l in [37j). We have checked 
this statement by high-precision calculations, which allow for "numerical proofs" in the 
manner explained in [HE! • The second term of ()A.14jl can be calculated analytically for some 
particular cases. For example, for f{n) = 1 we have 



^ /2r2\ x" 
^1 V y 



a+2-ra+l 



|) {l}a+l 

{2}„+i 



2" 



a-2 

E 



2 cos - 
2 



2(J + 1) 



ln^+2 

Lsj+2 (tt 



' (a-2)!^oV J 
ln^+22 



(\n^/x^ 



a-2-j 



X 



Lsj+2 (tt) (2cos^ 



(A.15) 



where we have used J2 ( a: ) T ^ ^ ^^(^ + x)? ^ 2, = 2 arccos , a = ±1 and Ls^ (6) 

k=l 
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is so-called log-sine integral [221 defined by 



■fc-i 



2 sin 



Ls, (9) = Lsf ^ (9) 



(A.16) 



The values of Lsj (vr) can be expressed in terms of (^-function for any j ^32jj. The sign a is 
defined via the relation x = e^"^^. The choice of the sign derives from the causal "+iO"- 
prescription for the propagator or as the sign of the square root of cos^ |. In the physical 
region of interest here we have a = —1. 

For X > 1/4, an non-zero imaginary part develops for the sum (jA.lBj) . It can be calculated 
in closed form: 



^ /2n\ 
Im y ^ 



n=l 



n y 



(a-2)\ 



j=0 



J 



J + 1 



Lsj+2 iTr-9)- Lsj+2 (tt) - 9 [ 2 cos 



On the mass-shell {x = 1,9 = = —1) we find the result 



9=2 arccos — i= 

V 43: 



(A. 17) 



E 



/2n\ 1 



1 ^ 



n J rr 



1 p 



)a— 1 



+i- 



a-1 ! 



{2}a+l 



In" 2 



(a -2)! a 



LSa ( - ) - LSa (tt) 



(A. 18) 



In general, however, analytical results for these types of sums are not expressible in terms 
of generalized Nielsen polylogarithms. For example, we have 



°° /9r)\ / 7r\ 00-2 ^2 arccos ^= 

2C2 - vr0 + J + In^ l2cos- 
2»-i g i^a-2^1 (lnv^^ri!2 ^ 



sm ^ 
d9 1 



cos ■ 



In \fx + In 2 cos 



9^ 



-\ a-2 



2012 {9) + 9\Yi V^coJ- 



fl2 

2C2 - + J 



ln-^(2cos^)+| 



J + 1 
Lsj+2 (tt 



Lsj+2 (tt) 



-2C2ln-'+^2 



Ls5i3 - ^) - Ls«3 (^) 



+ 



J + 1 
J + 3 



(2cos-j -ln^+=^2 



2i(T OI2 {9) ln^+^ (^2 cos - Lsc2,,+2(^^) 



-icr 

J + 2 



Lsj+3 (tt) - Lsj+3 (tt - 0) + ^ln^+^ ( 2 cos 



0=2 arccos —h= 

v4x 



X 



(A.19) 
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where J2 y ^ j = 2Li2 (x) + In (1 + x) and LsCjj (6') are new type of log-sine integrals 

introduced in |37j (where the properties of these functions are given in Appendix A. 2) and 
defined by 



LsCjj(6') = — / d01n 



2 sin — 
2 



In^- 



2 cos — 
2 



(A.20) 



For i,j > 2 these functions cannot be written in terms of Nielsen polylogarithms, but are 
related to harmonic polylogarithms which have been considered in [HHl • 

Using the results for the binomial series given above, we are able to calculate analytically 
the required first few coefficients of the e-expansion of hypergeometric functions of the 
following type 



I + biE, . . . , I + bjE, 1 + aiE, . . . , 1 + aK£, 2 + die . . . , 2 + (Ils 
2 + fie, ...,§ + fje, 1 + eiE, . . . , 1 + crE, 2 + ciE, ...,2 + cq-j-.rE 



qFp 3 r / ^ 3 



z] . (A.21) 



Rewriting this function as an infinite series and using the well-know representation 
T{j + oe) 



r(l + aE) 



U - 1)! exp 



-E 



k=l 



k 



SkU - 1] 



we obtain (for details we refer to Appendix B of ^7] ) 



{§ + f,E}i {1 + e,£}^, {2 + 



4z 



2z 



ttK+L~R 
■'-'■8=1 



(i + c,£)nii- 



1 + diE) 



(1 + 2brE 



^nL,(i + 2/..)f:(^')-^ 



X A , 
(A.22) 



where P = K + L + J and 



A = exp ^ 

\fc=i 

1 



k 



SkTk + 2 UkSk + 



1-eISiTi + 



'2f 



Y2 + 



Y^ 

3 



2UiSi 



SiTi + 2UiSi + 2UiSiSiTi + 2U2S2 + 2UlSl 



+ \s2T2 + \sItA + o{e^) 



(A.23) 



Here, we introduced new constants Ak, Bk, Ck, Dk, Ek, Fk,Tk,Uk: Ak = Z^ili^f, Bk = 
E.t/&^ Ck ^ ES^-^-'c^ D, ^ Etid^, E, ^ Et.et F, ^ ELift, T, ^ B, + 
Ck + Ek — Ak — -Dfc — -Ffc, Uk = Fk — Bk, = Ck — D^. With the expansion just worked 
out, analytical results are available for the first three coefficients of the e-expansion of the 
hypergeometric function ()A.21|) . In particular, we find 
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^ M 2-3e 



1 l-3el + x 



+2e'\ 3 Li2 (x) + 2Li2 (-x) + ln(l - x) + 2 ln(l + x) 



Y^-2|ln(l-x) + 21n(l + x)|e 



(A.24) 



For the analytical representation of the ^-expansion of V-type integrals we need also 
sums with shifted arguments (see Eqs. ()3.13p and ()3.16p ). The following two types of sums 
have to be considered: 



ai,...,ap; bi,...,bq;c\ ) 



jll,...,lp; Jl,...,Jq / \ 

ai,...,ap; fei,...,feq;c\*^/ 



= E C^) ^[SaAn-ir ■ ■ ■ [Sap{n-1)T^ [SbA^n-ir ■ ■ ■ [^^(2n-l)] 



n=l 



n / 



E 

n=0 



X 



n / [n 



^^[SaAnT ■ ■ ■ [Sap{n)X^ [SbA'^n)r ■ ■ ■ [^6,(2r^)]^' . (A. 



Cases when there are no sums of the type Sa{n — 1) or Sb{2n — 1) on the r.h.s. of ()A.25j) we 
characterize by a "—"-sign replacing indices (a,i) or {b,j) of /, J, respectively. 
With this notation we may write the following relations: 



2x 
1 



-^a;l;c(^) — 2^a;l;c(^) '^a;l;c+l (^) + ^1- 



Tr>^ 



2x 
1 

2x 



;c+l 



, [-^-;2;c(^) -^-;l;c(^) — 2</_;'2;c(^) 2J_.';^.^(x) ^-;'2;c+l(^) + J-;l;c+l{^) '^J--l;c+l} 



E L t:^^-^"" - 1) = E I I ^" 



2x \ n J n' 



n=l 



n 



-;2 

2KH 

(n + 1)^ (n + 1)^+1 



K(n) 



(A.26) 



B Fermion corrections to the pole masses at one-loop 



In this Appendix we present, for completeness, the well know one- loop relations between 

pole- and MS-masses of the gauge-bosons. We divide all corrections into bosonic (diagrams 

without any fermions) and fermionic (diagrams exhibiting a fermion loop) one's: Xy"* = 

X*^^] + X^^\ ■ where the purely bosonic contributions are given in Appendix B of 
y, boson yjermion f j o 

I. Using the notation 



V 



ml, 



1 + 



X 



(1) 



IGvr^ sin^ 6w , 

we may write the fermion corrections in the following form^^ 



(B.l) 
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For simplicity we assume a diagonal Cabibbo-Kobayashi-Maskawa matrix. 
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where 



Bo{m^,m^;mJ 




Bo{ml,ml;n 



quark ''^w''^H 



lepton ''^W''^H 



dxhii —^x H ^(1 - x) -xil - x) - iO 



denotes a scalar two-point function. 
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C Mass renormalization contributions 



For the parameter renormalization of a two-loop amplitude one requires the first order 
derivatives with respect to all relevant parameters of the one-loop amplitude as may be seen 
in ()4.4(i|l . Since, we are interested in the on-shell amplitudes mass derivations implicitly 
involve a differentiation with respect to the external momentum, as the "on-shell momen- 
tum" has been given the value of a mass. We restrict ourselves to consider the effect of 
the mass renormalization contribution, since for the charge renormalization the derivative 
with respect to the latter is trivial and it is included in ()4.46|) as a separate term. We thus 
consider 

i 



Di 



' dm'j 



(C.l) 



where the coefficients 6j are the one-loop mass renormalization counter-terms, which depend 
on the renormalization scheme. For the MS scheme they have to be identified with the Z^^2^^ 

of ()4.46|) . The explicit expressions presented below are written down for the third fermion 
family (z/^, r, t, b) in the approximation of vanishing r- and 6-mass. Accordingly in ()C1|1 . the 
relevant masses are indexed by j = W, Z,H,t. For simplicity we give the results assuming the 
Cabibbo-Kobayashi-Maskawa matrix to be the unit matrix. The contribution of a massless 
family may be obtained by putting rrit = 0{Ao{mt) = 0). The results read 
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D Bare two— loop contribution of the massless fermions 

As one of our results we present the exact analytic two-loop contribution to the on-shell 
self-energies of the gauge bosons. After reduction of the set of basic integrals to a minimal set 
of master-integrals by means of Tarasov's recurrence relations ^1] we obtain the following 
expressions: 
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with c = mw/mz — cos Qw and P^'^^ (p, m) = l/(p^ + m'^y and set {A, B, C, D, E} denotes 
the particles. The external momentum belongs to the mass shell. 
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E MS vs. pole masses at two-loop order 



After expansion of the diagrams with respect to sin^ 6w we get rid of one of the boson masses 
and write the functions Xy^ introduced in (jS.ip in the form 
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In particular for the Z boson propagator we ehminate mw and vice versa. Consequently, 
the coefficients in the above formula are functions of the Higgs and top masses and one 
of the boson masses. We expand this function with respect to my/m?m 

j=0 

and calculate analytically the first sixth coefficients This is not a naive Taylor expansion. 
The general rules for asymptotic expansions |28j allow us to extract also logarithmic depen- 
dences, or in other words, to preserve all analytical properties of the original diagrams. In 
the result of the asymptotic expansion all propagator diagrams are reduced to single scale 
massive diagrams (including the two-loop bubbles). In the finite part, we meet the following 
constants and functions: 
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where Nm is a number of massless fermion families, and ^{z) is the finite part of two-loop 
bubble master integrals defined in ISIl- We rewrite is as follow (see p?H]-P7j): 
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In parameterization HE.1|) _Bj = for j > 



43 



Below, we present the coefficients for Re A^jermion including fermionic contribution only, 
which are defined as 

ijull i,boson ~ i,fermion 

and AYhoson given in Appendix D of |1J. In the following we present the coefficients in the 
form AYj = BYj{m]j,mj) {m^y . 
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+9 In(cut) \n{ml^)tH'' - 72 In(cUi) ln(m^)tH' + 27^^' - — + ^^h^ - i^F (E.l) 

32 64 32 



177 339 

= l/Ah2UJHtH'^ + — ln(m^)tj/^t^^^ - —\Yi{m^^)uJhtH - ^FuntH^ 
2 16 

925 

\n{uJt)^HtH - 3/41n(cjH) ln(cjf)c^H - 2/3 \n{ujH)i^HNm + 3/41n(c^t) \n{m^)ujH 



48 

113 2 

+— (in(m^)) uJhtH + 16/3cj//tH^A^m - l/2Li ^(u;^)^^//^^ - l/4Li ln(u;H)^/f%^ 

37 37 
+— Li \n{ujH)ujHtH^ + I/2L1 \n{uJt)ujHtH + l/4Li ln(t^t)t^j^t^^^ - — Li ln(u;t)u;^^tH^ 

79 11 11 977 

+3/8L1 \n{LJh)LJt - — \Yi{ujH)uJHtH + Y^Li ln(c^t)c^H - ^TT^^H + l/167r^a;t + J^^htn 

1451 .27^. 1621 , , , 2 11 w 2 ^ , « ,w ,2 



-u;//tH - -ujHNm — M^t)^HtH + ln(m,y )a;/^ A^^ + 6 (ln(u;t)) uJntH' 



24 9 24 ' ' 6 

33 o 11t , / N 41 ^ 2 23 43 / , 2 x\2 2 
-zr^i^HtHSo - — Li ln(u;;i)u;H + zTzFuntH + 777:^^// - 77 (In(m^y) untn 
16 12 12 36 2 ^ ^ 

/ \2 7 99 

+ ln{ujH) \n{m^)ujH - (in(m^) j Uh - —Tc^uJhtH - 2/3ujHtHNm + —uJhtH^So 
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277 5 
-3/8Li ln{ut)ujt + — ln{ujH)u!HtH'^ + 3/16 {ln{ujt)f untn - —Fcvhtn + l/2h2UJHtH 

-^L2UJHtH^ + ^^{uh) ]n{m'^)u;HN^ - (ln{m'^)y uhN^ - l/24:ln{uH)ujH 

-l/2\ii{u!t)i^H + 5/Q\n{m'^)ujH - 5/6Fuh + 3/16Fut + ]^L2UJh - 3/8L2Ut 

-5/4\ii{ujh) ln{uJt)uJhtH + 5hi{ujH) ln{uJt)LUHtH'^ + 5/41n(a;H) ln(m^)a;j^t// 

99 99 
-5 ln{ujh) ln{m'^)ujHtH'^ - — ln{ut)ujHtHSo - 4:/3ln{uJt)ujHtHNm + — ]n{uJt)ujHtH'^So 

o Z 

29 31 
+16/3ln{uJt)'^HtH^N^ - — \n{LJt) \n{m^)uHtH + — ln(a;t) ln(m^)a;//t//^ 

99 99 
+— ln(m^)a;iftif5'o + 7/31n(m^)a;i/tijA^^ - — ln(m^)a;ijiif^5'o 

— — ln{mlir)uJhtH'^Nm - 2 (in(m^)) UhtRNm + 8 (in(m^)) iOhtH^Nm 

+2 \n{ujt) \n{ml.)ujHtHNm - 8 Hut) ln(m^)u;j^i^^2 (E_2) 



^S^2 = -^^LOni-l'l^mUuJHtH^ - 95904 ln(m^)a;HtH' - 98172 ln(m^)a;j^tj/ 
lU3o8 y 

+139968Fu;jfijf^ + 162u;j/y - 5292 ln(u;t)u;^/t// - 17280 ln(u;ft) hi{uJt)uJH 
-26880 ln(a;jj)a;jjA^^ - 97632 ln(u;t) ln(m^)a;ff - 139968 {hi{ml^)f untR 
-6480Li \a{ujH)uJHtH + 12960Li \ii{ujH)uJHtH'^ - 28080Li \ii{ujH)oJHtH^ 
+6480Li \{i{uJt)ujHtH - 12960Li hi{uJt)ujHtH^ + 28080Li \n{uJt)ujHtH^ 
+1296Li \{i{ujH)uJt + 39312 hi{uJh)ujHtH + 4320Li \ii{uJt)ujH - 120247r^a;jj 
+2167rV + 276426a;/fiH - 516672ujHtH^ + 1144080a;jfA^^ - 866052a;jf,52 
-1598A\n{u;t)u;HtH^ - 666816 ln(m^)a;i/A^^ - 4320Li \n{u;h)u;H - 69552Fu;HtH^ 
+886115u;if + 3888a;t + 559872 (in(m^))^ - 116424u;if,So 

+23040 ln(o;j^) \n{ml,)ujH + 133056 (in(m^))^ cuj^ + 155527r'^ untH - 1296Li ln(a;t)wt 
-2U7(}4\n{uH)ooHtH^ + 69984 {\n{ut)f untn + 276A8FuHtH + 6480L2CUjytH 
+28080L2Cc;//t//-'^ + 23040 In(tUi^) ln(m^)cUi^A/'^ + 114624 (in(m^))^ cUiyiV„ 
-20742 ln(u;H)w// + 287154 ln(cjt)cjjy - 608340 ln(m^)cjj^ - 7344Fcjh + 648Fo;t 
+432OL2CU/, - 1296L2CJt + 34992 \n{ujH) ln{u!t)uhtH - 93312 \n{ujH) IniutpHtn'^ 
-34992 ln{ujH) ln{mly)ujHtH + 93312 In(cjH) ln(m^)cuHtH^ + 69984 ln(a;t) ln{m^)LJHtH 
-559872 In(cut) ln(m^)tUiytij2 - 55296u;i/53 + 43008u;,,Ar„2 In(m^) 

-18432u;jjA^^2 (in(m^))^ - 6156^;^/ ln(u;t)y + 103680u;i/C(3) - 1296u;j^ {H^h)^ 

+23760u;jf 0.n{u}t)f - 20A80uhNJ^ + 29161n(u;H)^^t - 29161n(m^)u;t 
+228096a;//A^^,5o In(m^) + 27648^;/^ ln(a;t) ln(m^)A^^ - 99792u;hln{u;t)So 
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-32256u Hln{ut)N^ - 42624u;i/A^^7r^ - 5184^;;^ 0.n{u}h)f + 6156^;^/ ln(a;jj)y 
+162uhH^H)y^ + 18432a;jjA^„V^ + 414720u;jjA^^C(3) - 162ujHH^t)y^ 

+228096ujhSo In(m^) - 38016u;//iV„S'o - 221184a;//iV^S'3 - 15396A8(JhNmS2^ (E.3) 

= —tHUJH (-4:7671^ - 97250 + 8259 - 9504^^5*0 - 3552tj^ + 2376 Infect )S'o 
576 \ 

+3361n(a;t) - 95041n(a;t)tjj,5o - 3552 ln(a;t)tj^ + 612 On{^t)f - 2304 {]n{uJt)ftH 
-5761n(u;t) In(m^) + 20161n(u;t) ln{m'^)tH - 23761n(m^),So - 6601ii(m^) 

+95041n(m^)tH5o + 9601n(m^)tH - 36 (in(m^))^ + 288 (ln(m^))%i^ j (E.4) 



^1^2 = - Y^'^^' (^-11287r2iV^ + 576 ln{ml,)N^ - 1296 (ln(m2,))' iV„ + 10368C(3)iV^ 

+5832 ln(a;t) ln(m^)ijj - 24300 \n{u}t)tH - 31104 ln(a;t)ijj2 + 5832 {\n{u}t)ftH 
+15552 lii(m^)i/^ + 77761n(m^)iif^ + 46656 (ln(m^))%/f^ + 32238^/^ + 5364,5o 
+1452A^^ - 46656 In(cut) ln(m^)tH^ + 39691n(cut) + 11161n(m^) - 8997r^ - 27216^^^ 

+ 1458^2 + 2304^3 + 1728C(3) - 342 (In(cjf))^ - 1296 (in(m^))^ - 11664 (ln(m^))%H 
+9721n(a;t) In(m^) - 63365oiV^ - 17496^2^^^ + 9216^3^'^ - 8029 

-4752 ln(m^)5o - 4752 ln(m^)5oA^m + 2700 ln(u;t)S^ (E.5) 

A^i = — —tH^^H I 4767r^ + 1872,^0 - 8619 + 8640t^^5o - 6528tij - 21601n(u;t)5o 

576 Y 

+2184 In(cJt) + 8640 \n{ut)tHSo - 6528 \n{ut)tH - 612 (In(cjf))^ + 2304 (In(cjf))^ 
+576 \n{ut) In(m^) - 2016 In(cJt) ln(m^)t^f + 2160 ln(m^)So - 1860 In(m^) 

-86401n(m^)ijjS'oXst/ + 91201n(m^)ij^ + 36 (in(m^))^ - 288 (in(m^))^ J (E.6) 



^J^2 = (^-17496 ln(u;0 ln(m^)i^j + 35641n(a;0*i^ + 139968 ln(a;t)ij^2 

-17496 {\n{uJt)ftH + 22680 ln(m^)tH - 69984 ln(m^)tH^ - 139968 (ln(m^^))^ t^j^ 
-70686t// - 63114^0 - 13708iV„ + 139968 \n{cut) ln(m^)t//2 - 9561 In(cut) 
-10632 In(m^) + 499171^ + 104976i^^2 - 113238^2 + 3888C(3) + 90 {hi{^t)f 
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+4752 (in(m^))^ + 34992 (ln(m^))%if - 74881n(a;t)5o - 19441n(a;t) In(m^) 
+25344 ln(m^)5'o - 1031525'oiV^ - imUSiN^ + 4752 (in(m^))^ + 5768%'^ N^^i 
+25344 ln(m^),SoAU + 1728C(3) A^^ - 9552 In(m^) A^^ - 228^ (E.7) 

E.2 The expansion coefficients for the Z 

g 9 g 5 

\o = - g M^t) M^h)th + g In(^ft) ln(m|)4 - - ln{zt) \n{m\)th - ^q^^^^^ 

75 57 105 

+l/27rV -—Hzh)th + ^Hzh)th^ + — \n{ml)th + — ln(m|)ift' 
ib io io 

-^ln(m|)t,3 - 361n(m|)t,^ - (ln(m|))%,2 - 9/2 (ln(m|))%,3 

+36 (ln(m|))' + ^ In(zt)^/. - ^ ln(^t)t,,' + y ln(^O^ft' + 361n(zt)^?^^ 

+ \ {Hzt)f t,-^ (ln(z,))2 _ ^ gg (i^(^^))2 _ 9^^^ 

8 Id 8 

135^ . 27^ o 81 117, ,,,,,, , 
" ~ Y " ~ 16^'^ ^ "8" ^^^'^ Hm%)th^ 

9 63 

+9 ln(^,) ln(m|)ife3 _ 73 ln(^,) ln{ml)h' + -i^^ ^ _^^3 (^.8) 



= -3/81n(2;t) ln(2;ft)2;^ + ^th^hr - l/2FthZh - ^ln{zh)thZh + 10 ln{zh)th'^Zh 

99 

+5/81n(2;,j)t,j2;,jr^ + 5/4:ln{zh) \n{ml)thZh - 5\n{zh) \n{m\)thZh + — \n{m^z)ihZhSo 

o 

579 121 7 7 49 

+A/Z\n{m\)thZhNm + -^^hZh ^th^Zh - -ZhN^ + — Mzh)zh + — \n{ml)zh 

ln(2;t)z/i + + + In(^ft) hi{m%)zh - (ln(m|)) + 'ijMHhZh 

33 99 / n2 

-— iftZftS'o + —thZhSo + ln(2;ft) ln(m|)2;ftA^^ - (ln(m|) j zuN^ + 16/3ift^2;hA^„ 

55 281 

-—hi{zh)thZhr - 2/'i\n{zu)zhNm - — hi{m%)thZh + 81hi{m\)thZh 
48 16 

111 / 2^ .r 113/ , 2n\2 43 / , 2 x\2 9 205, , , 

+— ln{mz)ZhNm + — [ln{mz) ) hzh - — (^ln(m^) j tft + — \vL{zt)thZh 

135 , , , . 2 I Q /o 1 M /' 2 \ I 1^ /I / \\2 . I Q (-1 ^^2 + 2 



2 ■ln(2;t)ift + 3/81n(2;t)ln(m^)z/j + — (ln(2;t)) + 3 (111(2;*)) i/^ 2;^ 

99 40 / \2 

-5Fth^Zh - — \n{ml)th^ ZhSo - y hi{m\)thZhNm - 2 (^ln(m|) j UzhNm 

+8 (ln(m|)) thZhNm + ^ ln(^t)i,j2;ftr - 5/81n(2;t)ife2;,,r^ - y ln(2;t)ife2;feS'o 
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99 

-A/?,hi{zt)thZhNm + — Hzt)th^ZhSo + lG/3ln{zt)th''zhN^ - 5/4:ln{zt) Hzh)thZh 

3T 115 

+5 \n{zt) ln{zh)th^Zh - 8 \n{zt) ln{ml)thZh + — \T^{zt) \rL{m\)thZh - ^^^hZhr 

+ ^FthZhr^ + 2 ln{zt) ln{ml)thZhN^ - 8 In(zt) ln{ml)thhhN^ (E.9) 



K2 = -Y^'z^Nj + ll/3z^SoN^ + ^Fz^r' + |^ Hz,)zn'r' 

112 5673 
-— \n{ml)zh''Nm^ - 3/27r\zh'' + 1/2 {Hzh)f z^^N^ + Hzh)zh''r 

27,, , ,,2 2 1865 ^ 2 25979, , , 2 297 20 
-^iH^t)) thZh -^^FthZh + ^^^\n[zt)thZh + —ZhS2N^ 

.on / N 2 2 121 2 2 487163 2 377 2 2 7729 2 2557 2 2 
+181n(2;/j)t?j^2;ft^ H n^Zh hzh H h z^^ H Zh r 

-22 ln(m|)^ft25o + y^3^ft'iVm + '^Tchh^Nm - iQC{^)zh^Nm - Hzh)thZh^ 



384 ' ^ 32 ' ^ 27 ' ' 9 

^ 1 / \ 2343 -1 / 2 \ 2 1 / 2 \ 2 2 34 T3 , o \ 2^7 

ln(2;ft)2;t2;ft + — — h\{mz)thZh + -r mlm^jift zu + -r— In^m^j^ft A^^ 



160 ' ' 160 ^ " " 4 ^ 54 

+ y (ln(m|))' W - 54 (ln(m|))%ftW - ^ (ln(m|))' V^y^ 

+y (ln(m|))^ ZhNm^ + ^ \n{m\)ztZh + y 111(2;^) \n{m\)thZh 

20 

-9 111(2;,,) In(m|)tft2;2ft2 - — 111(2;;,) ln(m|)2;ft'Ar^ - 22 ln(m|)2;ft'5oAr^ 

y 

27 27 

— — hi{zt) \a{zh)thZh + 91n(2;t) hi{zh)thZh - — ln(2;t) hi{m\)thZh 

+541n(;2i) ln(m|)U'^ft' - 4/31n(^i) \n{ml)zh^Nm + 16/3^3^^' - '^^^^^^h 
49 2 3 297 2^ 55 2 o 7969 2,. 160 2,. 2 ^./on 2 

+ Y^^ftV=^ + -^^/^ -^2 + -^Zh^So - -T^Zh^Nm + —Zh^Nj - 7C(3)^ft' 

8413, , , 2 5 „ , ,,2 2 50209, , 2 n 2 77/, , 2n\2 2 
+ ^^ln(^,,)^,, + — {\n{zh)) Zh +-^ln(m2)^ft - — (In(m^)J Zu 

38411, . . 2 35 , ,,2 2 3 , ^ 2 17,, 2 2 



1440 ^"(^*)^'^ + (^"(^*)) ^'^ - i^^*^'^ + V + ^ In(zi)^/.'^/.' 

^^'^^ ■ ln{zt)zh^r - ln{zt)zh^r^ + ^ ln(^t)^,i^r^ - ^ ln(2;t)^ft^r^ 



10240 ' ^ 23040 ' ^ 384 ^ " 32 

55 14 113 

+ — \n{zt)zh^So + — \n{zt)zh'^Nm + 5/6 ln(2;i) \n{zh)zh'^ + — \n{zt) \\i{m\)zh 

53695 ^22 27 ^ , 2 89 ^ 2 1441 ^ 2 2 241 2 3 x 
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165 

Af^ = ln{zt) \n{zh)zh + l/2\n{zh)thZh —ln{ml)thZhSo - 8/31n(m|)tfe^/iiV„ 

9919 607 14 

-thZh + —thZh + —ZhNm - l/3\n{zh)zh - 2 \n{ml)zh + 7/6ln{zt)zh 



288 3 " " g 

+l/4Fi;,2^fc _ 2/9^^ _ 2 in(z^) ln(m|)^ft + 2 (ln(m|))%ft - 3/87r\zh + '^hZhSo 

2 32 

-Wbth^ZhSo - 2 In(zft) \n{m\)zhNm + 2 (ln(m|)) z/jiV„ - —th^ZhN^n 

12()7 790 
+31n(2;,j)ife2;fer + 4:/3]n{zh)ZhN^ + ]n{ml)thZh — \n{m\)thZh 

127 2 2 625 

-ll/'i\Yi{m\)zhN^ — (ln(m|)) thZh + 57 (ln(m|)) t/,^^;/, - — \n{zt)thZh 



8 ^ ' V ^ 12 

+— ln(2;t)ih^2;ft - ln(2;t) \a{m%)zh + — (ln(2;t))^ t^z;, - 12 {hi{zt)f thZu 



80 2 
+165 ln(m|)t,,2^,,5o + y hi{ml)th^ZhN^ + 4 (ln(m|)) t,,2;,,A^^ 

2 165 
-16 (ln(m|)) thZhNm - 3 ln(2;t)t^2;ftr + — \{i{zt)thZhSQ + 8/3 ln(2;t)i;j2;,,A^^ 

32 29 
-l<Qbhi{zt)th zuSq - — hi{zt)th ZhNm + — hi{zt) hi{m\)thZh - 451n(2;t) hi{m\)thZh 

+3/AFthZhr - 4 In(^t) \n{ml)thZhNm + 16 In(zt) ln(m|)ife2^ftAr^ (E.ll) 



Af, = 16/37r2^,27V^2 + ^z^SoN^r, - ^ ln(z,)^.V2 + f ln(m|)^,27V„2 

9 54 o 

1321 391 

+37i\zh^ - {ln{zh)f Zh^Nm - T7— Hzh)zh^r + 9 {Hzt)f hzh^ + —FthZh" 

iUoU zoo 

17969, , , 2 „ ^ 2^ 167 2 2 1913 n , 22 233 2 

ln(zt)i^z^' - 720zfe'>S27V^ - -^T^ - -^^hZh - ^'^th^zn^ - ^^ftV 

+ + ^ HmDzn'So - ^^S^Nr. - ^^t^'z^N^ + 176C(3);^,^iV^ 

ln(2;?,)i^2;^^ + y ln(2;fe)2;,jV^ - ^ ln(2;?j)2;^^A^^ + y ln(2;^) \n{m\)zh 

509 595 
+l/51n(^/,)^t2;^ - — ln(m|)ift^ft^ + 12\n{m\)th zh — \n{m\)zhNm 

-18 (ln(m|))%,z,^ + 72 (ln(m|))%,^z,^ + ^ (ln(m|))%,2iV^ 

-16/3 (ln(m|))' Zh^N^^ - 1/5 ln(m|)^t^ft + ^ ln(^^) ln(m|)^^2^^ 
242 

\ri{m\)zhSQN,n + 9 ln(2;t) \n{m\)thZh' - 72 In(zt) ln(m|)t?j^2;/i^ 
40, , xw 2x 2.r 208^ 2 475837 2 -,0^ 20 289 3^ 
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-r^z^N^ - ^^z.'Nj + ^^mz,' - 1^ Hzu)z^ - 1/2 (ln(z,))' 
ln(m|).,^ + 41 (ln(m|))' + ^ ln(.,)^.^ - 7/6 (ln(.,))' 

15 1U8U 54 

46, , , 2 3 619, , , 2o 116, , , 20 ,,,,,, 2 

— ^ \n[zt)zh f ^'^{.zt)zh ^0- ^ i-Azt)Zh Nm- — hi{zt) ln{zh)Zh 

Hz,) Hml)z,' + ^^Ft,W - ^Fzn'r - ^Fz^r' + (E.12) 



36 ' ' ^' " 72 " " 288 12 18 



183 32 

= -4/31n(zt) \i\{zh)zh - Ij'iHzi^thZh + — ln(m|)Uz/iS'o + y ln(m|)U2;/iiV„ 

1567 2138 2 56 . , , / x / , / 2 x 14 , , x 

g"^'' ~ Ori ^ 4/91n(2;^)2;^ + 8/3111(7^2)2;^ - — \n[zt)zh 

-l/?,FthZh + ^^ft + 8/3 \a{zh) ln{ml)zh - 8/3 (ln(m|)) 2;^ - ^th2;/j5'o 

2 128 

+183t^'zfe>5o + 8/31n(zfe) \n(ml)zhN^ - 8/3 (ln(m|)) ^^.A^^ + —th^ZhNm 

/ X 16, . X 518, , 2x 2324, , ox 2 

-41n(2;,,)t/,2;,,r - y \n{zh)ZhN^ — ln{mz)thZh + m{mz)th Zh 

44, , 2x 31 /, , 2n\2 62/ , 2n\2 2 524, , , 

+ "9" + Y (Inlm^) j t/,^/, - y (^111(771^) j t/j + — ln(2;t)i;j2;,, 

2042 , , , . 2 I /I /Q 1 / \ 1 / 2 x c /Q /'I /' ^^2 + 1 ^2 . ^ ^2 , 2, 



ln(^t)ift + 4/3 \n{zt) ln{m'z)zh - 8/3 (In(^O) ^/.^/^ + ^ (In(^t))' h^Zh 



-lS3]n{ml)th^ZhSo - ^ln(ml)th''zhNm - 16/3 (ln(777|))%^z^A^^ 

64 2 183 32 
+ y (ln(77i|)) thZhNm + Ah\{zt)thZhr — h\{zt)thZhSQ - — \n{zt)thZhNm 

128 

+183 ln{zt)th^ZhSo + — \n{zt)th'^ZhN^ - 5/2 \n{zt) \n{ml)thZh + 10 111(2;^) \n{ml)thhh 

64 

-Fi;,z^r + 16/3 ln{zt) ln{ml)thZhNm - — Hzt) ln{ml)th^ZhN^ (E.13) 



-3/2n%zn' + 4/3 (ln(;.,))2 z.^iV^ + ^ ln(z,);.,V - 9/2 (Inizt))' t,z^ - ^F^zh' 

oiU Zib 

34453 2^ 2655 2^ ^ 995 , 2 ^ 4208 2 ^ ^q, 2 2^ 233 2 
■ \n{zt)thZh H r— -z/i S2Nm + — TT + -—thZh + 39th Zh + -j^Zh r 



540 V ^/ " " ■ 2 " ^ 216 135 » » ■ 

92 22 638 2 x 2^^ 1408^ 2., ^ 1025 , , 860 2., 
oi y y o4 o 

577 , . 2 184, / x 2 s 560 2 at 160, , . 2\ 2 

ln(2;^)i,j2;/j - — \n[zh)zh ^ + ln(2;ft)2;ft - y^ 111(2;^) ln(77i2)2;/j 

50 



+9(ln(m|)) W-36(ln(m|)) V - y (ln(m|)) ^^^TV^ 

2 4 160 

■16/3 (ln(m|)) zhNJ^ + — ln(m|)ztZ/i - — hii^Zh) ln{ml)zh'^Nm 

)38 

— ln(m|)^^'5'oA^^ - 9/2 In(zt) ln(m|)ifeZ;,' + 36 ln{zt) ln( 
L12, , , 2n 2.r 352^ 2 1126103 ^ 2655 



638 

-— \n{ml)zh'SoNm - 9/21n(zt) ln(m|)ifeZ;,' + 361n(zt) ln(m|)i^'^^2 
112, . 2^ 2., ^ 352 2 1126103 2 _^ 2655 2^ 1744 2c, 

68035 2,, 1040 2,, 2 533 2 ^ 2192 2 ^^fonf 2 

3458 2 ^ 2 113/, . 2^^2 2 206869 2 ^ 23 2 2 

—ZtZh + 2/2,Fzh^ + 601n(^0^/.''^^' - ^ H^W^ - H^t)zh^r'' 
184, , , 2 3 983, , , 20 64, , , 2.r §0 ,,,,,, 2 

ln{z,) HmDzn' - ^-^Ft,W + ^Fz^r + \Fz^r' - f^Fz^r' (E.14) 
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Figure 4: The dependence on the number of coefficients of the expansion with respect to 



sin^ 6w of the two-loop corrections S = I H^'* + I!^''!!^''' \ (see 4-4(^ ^■^ o, function of the 



(i)tt(i)/ 



MS 



Higgs mass. The dotted, dashed, dot-dashed and full lines show results obtained with the first 
one, two, three and all calculated (six) coefficients, respectively. Upper plot: for intermediate 
Higgs masses. Lower plot: for heavy Higgs masses. 
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Figure 5: The dependence on the number of coefficients of the expansion with respect to 



m\/mf (i = H,t) used for the evaluation of the two-loop corrections. We show 6 
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and full lines show results obtained with the first one, two, three and all calculated (six) 
coefficients, respectively. Upper plot: for intermediate Higgs masses. Lower plot: for heavy 
Higgs masses. 
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Figure 9: Corrections to the relation Az = Mz/m\{Mz) — 1 as a function of the Higgs mass 
Mh for heavy Higgs masses. Upper plot: the various two-loop corrections. Lower plot: the 
complete one- and two-loop correction. 
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Figure 10: Two-loop imaginary part 6 = Imsp^y/my (imag) and absolute value 6 = 
\sp,v\/'>TT'v (abs) of the pole position sp for the Z (Z) and the W (W) as a function of 
the Higgs mass Mh- Upper plot: for intermediate Higgs masses. Lower plot: for heavy 
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